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Abstract. We investigate a possibility of heating of the loops and other closed magnetic structures in active regions 
of the solar corona by the flow of solar wind (plus other flows that may be present) across the magnetic field 
lines (that are perpendicular to the flow), in a similar manner as a conventional MHD generator works. A simple 
model is formulated which allows to calculate the typical currents generated by this mechanism. This enables to 
make a judgement whether the loops in solar active regions can be heated by the proposed mechanism through 
the associated Ohmic dissipation. It is shown that even if the heating region width is of the order of a few proton 
Larmor radii (this effectively means that a loop is represented by nano-threads, by analogy with the nano-flare 
coronal heating mechanism, with width of a few proton Larmor radii) and the plasma flow speed is 50 km s“^ then 
only about 1% of the heating requirement of coronal loops could be met under this conjecture. Other implications 
of the model, such as a predicted spatial heat distribution along the loop as well as application of the model to 
the other places e.g. chromosphere and corotating interaction region structures in the solar wind are discussed. 
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1. Introduction 

The coronal heating problem, the puzzle of what main¬ 
tains the solar corona 200 times hotter than the photo¬ 
sphere, is one of th e main outstanding questions in so¬ 
lar physics (see e.g. lAschwandeiJ (1200^ for a recent re¬ 
view on the subiectl. lAschwandenI ll20o'3l categorises the¬ 
oretical models of coronal heating processes into 5 main 
groups, according to the main underlying or driving phys¬ 
ical processes. These groups are (i) DC (Direct Current) 
stressing and reconnection models, (ii) AC (Alternating 
Current) wave heating models, (iii) Acoustic heating, (iv) 
Chromospheric reconnection, (v) Velocity filtration. The 
main focus of the solar physics community research in the 
coronal heating theory has been concentrated on the DC 
and AC type models. 

These two types describe the electro-mechanic coronal 
response to the photospheric driver that provides the ul¬ 
timate energy source for heating. It is believed that the 
perturbations of magnetic field are generated e.g. by the 
convective motions in the solar convective zone. Then the 
perturbations propagate through the photosphere to the 
corona and ultimately constitute the main energy source 
of coronal heating. These perturbations propagate with 
the Alfven speed. The judgement whether the heating is 
of the DC or AC type can be made by comparing time- 
scales of the heating with the Alfven transit time along a 


coronal loop: (i) When the photospheric driver, e.g. ran¬ 
dom motion of magnetic field line footpoints, changes the 
boundary condition on time scales much longer than the 
Alfven transit time along a coronal loop, the loop can ad¬ 
just to the changing boundary condition in a quasi-static 
manner. Hence the coronal currents are almost direct ones, 
which attributes a model to the DC type, (ii) However, 
when the photospheric driver changes faster than a coro¬ 
nal loop can adjust (e.g., by damping and dissipation of 
incoming Alfven waves and fast and slow magnetosonic 
waves generated due their coupling to Alfven waves on 
plasma inhomogeneities), the coronal loop experiences an 
alternating current, which attributes a model to the AC 
type. 

The DC and AC coronal heating models are subdivided 
further by possibilities of how the cur rents are dissipated: 
eithe r by magnetic r econnection jSturrocl^^Uchi^ 


JPark^ 1983albL IPriest et al ] l2nn2L IPrieslI 1200, it 
1 l2nn■3^. Ohmic d i ssipation via current 
casca ding llBallegooiienI Il98fil. ICalsgaard fc Nordlundl 
199dl. and viscous tur bulence IIHevvaerts fc Priestlll992L 


19811 

Priest et al 0^ 


Einaudi et alJll996bla^ in the case of DC models, or by 


Davilalll987l IPoedts et al. 

ll989llOfma,n et a,lJll994 199,5. 

Erdelvi & Goossenslll994 

199.5111994 

VasQuez & Hollweg 

2004 

.iR.udermanl 1999^, phase mixing ( 

Hevvaerts & Priest 

1983 

. iNocera et al.lll98fil IParkeilll99 

1 iNakariakov et al. 
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19971 iR.uderman et al.lll998lll99t 

.iDeMoortel et al, 

200C 

Botha et alJl200nllTsiklauri et al. 

l200lllHood et al. 

2002 

Tsiklauri & Nakariakovll2002l iTsiklauri et al.ll2002l 

2003 

2005albr), and turbulence Jlnveraritv & PriestI 

1995 

Dmitruk et alJl200lh in the case of AC models. 


As an alternative to current dissipation, some heating 
could also be produced by compressional wav es (i.e., by 
acoustic waves or shocks llKunerus et al.lll9^ b It should 
be mentioned that a problem which models with com¬ 
pressional waves face is that such waves cannot propagate 
from photosphere into corona due to the density gradi¬ 
ent. However, interest in these have been recently renewed 
due a possibility of generation of fast and slow magne- 
tosonic waves vi a their coupling to Alfven waves on plasma 
inhom o geneities[ Nakariakov et al.l ll997|) , 
'20f^. jT^klaurrfcNakariaim ll2002ll . 

2nn2[l2nn,1) . 

In addition yet another approach is adopted by the 
chromospheric reconnection models: the idea is that 
chromospheric reconnection through generation of frag¬ 
mented flux-tubes, strong currents, fast magnetosonic 
waves, and upw ards plasma flow s can contribute to th e 
coronal heating I Sturrock I lll999ll. iR.vutova et alJ ll200l| l. 


E 


Tsiklauri et al 


Tsiklauri et al 


iirusawa 


fc SakaiU2000l).ISakai et alJ 


I2nnna 


I2nnir 


Also a completely different physical mechanism has 
been proposed named velocity filtration, which is based 
on the influence of the gravitational potential field in the 
corona on a postula ted non-Maxw ellian chromospheric ve¬ 
locity distribution I Scuddei)ll994li . 

In this paper we investigate yet another mechanism for 
the heating of solar coronal loops and other closed mag¬ 
netic structures in general. We study a possibility whether 
the loops could be heated by the flow of solar wind plasma 
(plus other flows that may be present) across them by 
generating currents in similar manner as a conventional 
Magnetohydrodynamic (MHD) generator works. In sec¬ 
tion 2 we formulate a simple model which allows to esti¬ 
mate typical currents generated by this mechanism. We 
show that if the heating region width is of the order of 
few proton Larmor radii and the plasma flow speed is 50 
km s“^ (which some observations do report) even then 
the heating of coronal loops could be only 1% of the heat¬ 
ing requirement under this mechanism. In section 3 we 
discuss some relevant aspects pertaining to the studied 
mechanism. 


2. Main Consideration 

Similarly to the stressed-induced cu rrent cascade models 
(described e.g. in i AschwandenI 1 20041) . p. 366), we try to es¬ 
timate energy budget necessary to heat the solar corona. It 
is assumed by the stressed-induced current cascade mod¬ 
els that the random footpoint motion stirs up a potential 
magnetic field, and hence non-potential fields and pertain¬ 
ing currents are generated. One can estimate a significance 
of the Ohmic (Joule) dissipation associated with these cur¬ 
rents and whether they can contribute to the coronal heat¬ 
ing in a following manner. The volumetric heating rate 






Fig. 1. The sketch of the model. 


Eh requirement for a loop with length I which is heated 
at both footpoints in the units of erg cm“^ s“^ is 


Eho = = 2 X 10-3 


I 


I 


10^°cm 


Here we use d Fnn 


= 10^ cm ^ 


-1 


,-i 


( 1 ) 


.. from Table 9.1 in 

lAschwandenI ll2004li l. The required current density J for 
the Ohmic dissipation, given by Eh = I is hence. 


^_ / I \-i/2 

J = ^/uEho = 1.1 X 10" 


116 g-l 


( 2 ) 


for a 


Here the classical conductivity of(T = 6xlO^ 

T = 2 MK corona has been assumed. 

In general the geometry of a magnetic field in active 
regions is quite complex. It mostly consists of closed struc¬ 
tures such a loops, arcades, canopies, etc. Instead of con¬ 
sidering curved magnetic fields, we simplify treatment by 
approximating active region loop with a straight slab of 
length I {I and associated with the loop magnetic field are 
both along x-axis, B — Bz) and radial thickness d {d is 
along x-axis). Note that d plays the same role as a length 
dimension along the pla sma flow) of a MHD generator 
iKrall fc TrivelDiec3ll973ll . The plasma flows along x-axis 
with velocity U = Ux. The sketch of the model is given 
in Fig. 1. For our model we can write Ohm’s law in the 
following form fe.g. lKrall fc Trivelniecel I WT^)): 


U X B 


= -E 


J 

a 


Ba 


-J X B. 


(3) 
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Here we use standard notation with E, J and a being 
electric field, current density and plasma conductivity re¬ 
spectively. Wee is an electron gyro-frequency and r is time 
between collisions. Note that we can omit Hall current, 
the last term, in the Eq.(3) because we assume that cir¬ 
cuit is open in the direction of plasma flow and no current 
can be sustained. However, current due to the Lorentz 
force can flow in y-direction (i.e. into the plane of Fig. 1. 
Note: the generated current flows not along the slab, but 
across it, in y-direction! This is quite different from the 
stressed-induced cu rrent cascade models (described e.g. in 
lAschwanderJ ll2004ll . p. 366) in which azimuthal twisting of 
the magnetic field lines of a cylindrical loop generates cur¬ 
rent along the loop. In our model current flows from one 
edge of the footpoint to the other edge of the same foot- 
point, but involving entire bulk of the slab (which mimics 
the loop)). In a steady state usual MHD equation of mo¬ 
tion 


dU JxB 

p-^ + pU -VU = -, 

ot c 

is reduced to the following form 

p dU^ _ \J X B\ 

2 dx c 


(4) 


(5) 


pC/2 


|J X B\ 


d. 


( 6 ) 


_ 1 pU'^c 
~ 2 Bd ' 


(7) 


Note that we neglected pressure gradient in the Eq.(4) 
because dominant force in the solar corona is associated 
with magnetic field, i.e. plasma /3 <C 1. We can simplify 
Eq.(5) by putting d/dx ~ 1/d 


Hence we can estimate the magnitude of the current den¬ 
sity J as 


Note that the same as Eq.(7) applies to the current density 
produced by and MHD generator. However, for the latter 
d is replaced by the length of the MHD generator along 
which a plasma flows. 

Based on Eq.(7) we can make following estimates for 
the current density . We use following typical values (e.g. 
lAschwanderJ J2004ll l of parameters: An electron number 
density of an active region loop = 2.0 x 10® cm“® (hence 
the density is p = pempUe- We put p, = 0.61 corresponding 
to the a totally ionised gas composed of 90 % of hydrogen 
and a 10 % of helium i Sarro et al.lll999l)'): In closed coro¬ 
nal magnetic structures. iNoci et alJ JIQQTII . using UVCS 
operating on board of SOHO, showed that the spatial dis¬ 
tribution of the r.m.s. velocities of OVI ions, in stream¬ 
ers, is characterised by the fact that in the brightest re¬ 
gions (the streamer core, around the axis, at the lowest 
heights) the r.m.s. velocity is the smallest. The r.m.s. ve¬ 
locity is roughly constant in the core while, outside it, it 
incre ases, both moving in radial distance and in heliolati¬ 
tude. iNoci et ^ jlQQTll measured r.m.s. speeds (in a core 
of a helmet streamer) to be of the order of < 50 km s“^ 


between 1.5 — 3i?Q. We use this value of 17 = 50 km s“^ in 
the Eq.(7), however this should be used with caution that: 
what is actually measured is the r.m.s. speed i.e. some 
bulk motions of plasma carrying OVI ions with it (through 
measurement of OVI A1037/A1032 ratio) which may not 
necessarily be the velocity of the bulk plasma (protons and 
electrons) across the magnetic field of the closed struc¬ 
tures. It should be also mentioned that so me observations 
(see e.g. bottom right panel of Fig. 4 in i Stracha.n et !ii] 
I2flfl2ll l indicate that outflow speeds of 0®+ ions for the 
two UVCS slit orientations at 2.33 Rq from the Sun’s 
centre i.e. 1.33 Rq above its surface, above a helmet 
streamer, are less than « 20 km s“^ in the ce ntral region 
of the streamer, which is in contra diction w i th iNoci et iii] 
(Ei^’s measurements. Although lOfmanI l|2000l) argues 
that outflow speeds for the bulk of the plasma (protons 
and electrons) should be higher than that for ions 
(because the latter drag out through Coulomb collisions), 
still plasma flows could be as low as ~ 10 km s“^ (Ofman, 
private communication). Adopting such low values of the 
plasma outflow speeds would diminish our estimates con¬ 
siderably as the current is proportional to the t/®. For 
our estimates, at any rate, we use more ’’optimistic” value 
of t/ = 50 km s“^ quoted above (also see our comment 
(ii) in the Discussion section). A magnetic field of 100 
Gauss. Note that we only use longitudinal magnetic field 
strength, i.e. we totally ignore azimuthal twist of the mag¬ 
netic field. There is only evidence for weak twists in coro¬ 
nal loops, of the order of B^/B^ « 0.01, from measure¬ 
ments as well as from the kink instability criterion. Thus, 
taking into account only B^ component of the magnetic 
field seems as a good approximation. This is quite dif¬ 
ferent from the stressed-induced current cascade models 
in which azimuthal twisting of the magnetic field lines 
of a cylindrical loop generates current along the loop. In 
our case the current flows from one edge of the footpoint 
to the other edge of the same footpoint, involving entire 
bulk of the slab (that mimics the coronal loop). An ad¬ 
ditional thought is needed when specifying the parameter 
d. Initially we can try p utting an observed loop width. 
lAschwauden et aTI ( 200(11 present data based on 41 active 
region loops. They found that typical loop width is few 
Mm. Hence putting loop radial thickness d ^ 1 Mm we 
obtain. 


J = 0.08 


IMm 


(esu). 


( 8 ) 


This estimate turns out to be about eight orders of mag¬ 
nitude less than the requirement to heat solar corona cf. 
Eq.(2). 

In the next estimate, following lAschwanden et aTI 
J 20 onli . we adopt the multi-tread model of a coronal loop. 
A good graphi cal representation of this idea can be found 
in Fig. 1 from lAschwanden et al ] tOOlt . in which a loop 
is represented by 100 fine loop threads with random foot¬ 
point positions and random fluxes. Thus assuming value 
for d — 0.01 Mm, our estimate of current density will still 
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fall six orders of magnitude less than the coronal heating 
requirement. 

It has been a general trend that as the resolution 
of EUV telescopes increased more finer and finer struc¬ 
tures were observed. It is a general agreement now 
that loops have much finer substructure than previously 
thought. The lowest possible scale for the MHD ap¬ 
proximation to be valid is few Larmor radii of protons. 
Which is tl = Vl/wcp. Wcp = eB/rUpC is a proton gyro- 
frequency. Putting V± = /nip with temperature 

of T = 2 MK and B = 100 Gauss, is then 23.23 
cm dShibata fc Tanuiii^ ll200lll give similar estimate for 

within a factor of -v/S)- Therefore assuming that loop 
consists of substructures, what we call nano-threads (by 
analogy with the nano-flare coronal heating mechanism)), 
with a width of a few proton Larmor radii for the gener¬ 
ated current we obtain 


J= 1.1 X 10® 



(esu). 


(9) 


This is only about 1% what is required for the heating 
of solar corona cf. Eq.(2), which obviously cannot be re¬ 
garded as a sizable contribution to the heating. By obtain¬ 
ing this, one should be also aware of other considerations: 

(i) One may question a possibility of existence of such 
nano-treads. There is no direct evidence in the observa¬ 
tions that the coronal loops consist of threads with width 
of a few Larmor radii. In fact, essentially all finest loops 
observed with TRACE have a width of order 1000-2000 
km and have a monolithic structure, that is a single tem¬ 
perature at any cross-section, which presently does not 
support the nano-thread concept. The TRACE CCD cam¬ 
era has 0.5 arcsec pixels which is 366 km on the sun. Hence 
it is possible that the smallest observed 1000 km wide 
monolithic structures are simply too close to the resolu¬ 
tion limit and future high spatial resolution space missions 
may or may not reveal further sub-structuring ultimately 
corroborating or disproving our nano-thread conjecture. 

Situation here is possibly similar to the stressed- 
induced current cascade model, which based on the ran¬ 
dom footpoint motion (by twisting the field lines) pre¬ 
dicts a current density, which is about 5 orders of magni¬ 
tude smaller than re quired to sati s fy the coronal heating 
requirement Eq.(2). lAschwa.ndenI (|2nn4li concludes then 
that the Joule dissipation is generally inefficient in the 
corona, unless 10® smaller transverse length scales can 
be produced. This was the main motivation of the model 
bv iBallegooiienI lll986h . who proposed a current cascade 
model, where free magnetic energy is transferred from 
large to small length scales in the corona as a result of 
the random motion of photospheric footpoints. 

(ii) The other mechanisms reviewed in the introduction 
section are main contributors to the heating budget. And 
MHD generator mechanism contributes only tiny 1% even 
at the smallest possible scale. 


3. Discussion 

The main novelty of this study is exploration of a con¬ 
jecture that loops in the solar corona could be heated by 
currents which are generated by solar wind (plus other 
flows that may be present) across the magnetic field lines 
in similar fashion as conventional MHD generator works. 
In previous sections we gave estimates of generated cur¬ 
rents and appropriate heating scales. 

We would like to close the paper by discussion of fol¬ 
lowing four points. 

(i) Spatial distribution of heat deposition along the 
loop holds the key to an understanding of the physi¬ 
cal precesses that are responsible for heating of solar 
corona in general. This quantity can be inferred e.g. 
from high resolution images produced by the Yohkoh 
Soft X-Ray Telescope, the SOHO Extreme-Ultraviolet 
Imaging Telescope a nd the Transition Region and Coronal 
Explorer ('TB ACEl. lAschwanden et alJ (|200nll established 
that long loops > 100 Mm are preferentially heated near 
footpoints with the heating scale height of about 20 Mm 
and are far from hydrostatic equilibrium. While shorter 
loops ~ 20 Mm are heated almost uniformly and are 
predominantly in hydrostatic equilibrium. Although the 
heating budget requirement is not met by our model, it 
predicts that loops with weak curvature (radius of cur¬ 
vature is much greater than loop length, i?curv ^ 0 
will be heated almost uniformly, as the most part of the 
loop (with its magnetic field) will be at 90° angle to the 
plasma flow. Thus, U x B term will be significant in all 
parts of the loop. This is actually what is seen in the 
lAschwanden et al~l ^l2000^ data: e.g. if we look at images 
of loops 2,3,4 (length ~ 20 Mm) of the TRACE 171-f 195 
data they all have large curvature radii (are almost flat) 
and are heated uniformly, i.e. in accord with our model 
prediction. However, as we estimated above heating re¬ 
quirement is still not met by our model. Hence, the dis¬ 
cussion about the spatial distribution of heat deposition 
along the loop has only academic interest. For the loops 
with strong curvature (radius of curvature is of the or¬ 
der of loop length, i?curv — 0 our model predicts that 
footpoint parts will not be heated as the magnetic field is 
almost parallel to the plasma flow and thus at footpoints 
U X B term will be small and only the apex part of the 
loop will be heated. If we look at images of loops 30,31,33 
(length ~ 200 Mm) of the TRACE 171-1-195 data they all 
have curvature radii comparable to the loop length (are 
quite curved) and are heated near footpoints. This faces an 
additional problem to our model which predicts no heat¬ 
ing at footpoints and heating concentrated at the apex 
for the reasons explained above. Thus, loops with strong 
curvature are probably heated by some other mechanism 
(near the footpoints). However, as for the coronal heating 
problem itself this does not pose a serious problem as the 
number of short (~ 20 Mm) and flat loops is more signifi¬ 
cant as compared the long (> 200 Mm) and curved loops, 
which are more rare in occurrence. 
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(ii) We also would like to comment on the value of den¬ 
sity used in our estimates of the magnitudes of currents 
generated by our mechanism. We used the observational 
value of Uf. = 2.0 x 10® cm“^. In a conventional MHD 
generator value of density is prescribed by the strength 
of the magnetic field across the flow (which tends to re¬ 
sist to the flow) and the strength of the gradient of the 
plasma flow through the generator. In the case of coronal 
loops however there are also up-flows of pl asma from the 
chrom osphere and transition region. Also, ISvestka, et a,1.l 
dlQQSll suggested that in active regions there are loop sys¬ 
tems that are embedded in hot coronal structures look¬ 
ing in soft X-rays like fans of coronal rays. These struc¬ 
tures ar e formed during the fl are and extend high into the 
corona. ISvestka et al.l lll998ll found that such structures 
create mass flow from the active region into interplane¬ 
tary space and these are the sources of a density enhance¬ 
ment in the solar wind. This suggestion is supported by 
synoptic maps of solar wind sources constructed from scin¬ 
tillation measurements which show a source of enhanced 
solar wind density at the position of the active region. It 
is not entirely clear at this stage what effect the above 
complex factors could have on our simple calculation, but 
we believe in the overall correctness of our model and its 
predictions. At any rate one can conjecture that the ob¬ 
served density in the active regions of rig = 2.0 x 10® cm“^ 
is caused by the fact that the closed magnetic field lines, 
transverse to the plasma flow, resist the flow making it 
dense, such that the observed density is obtained. The 
same occurs in an MHD generator. 

(iii) We would like to point out the difference between 
our model and the AC and DC models mentioned in the 
introduction. For the latter two the main source of energy 
is the convective motions that drive either reconnection 
events or generate waves which then dissipate. In our case, 
however, the source of energy is in the plasma flow across 
the magnetic field lines, i.e. solar wind plus other flows 
that may be present. Commenting on possibility of actual 
flow across the magnetic field in the coronal conditions we 
would like to stress that the parameter which controls the 
possibility is a, the plasma conductivity. Since we do not 
use any enhanced values for a and the bulk flows in th e 
cores of helmet streamers are observed llNoci et alJll997ll . 
we believe that the flows across the magnetic field are pos¬ 
sible. Note that smallness of plasma-/? in the corona is not 
a problem for our model as this parameter controls only 
compressibility of plasma, not the cross-field flow control. 

(iv) Also, we would like to point out a possibility to 
apply our mechanism in the chromospheric footpoints of 
the loops or in the corotating interaction region struc¬ 
tures in the solar wind. For example, in the chromosphere 
the neutrals are abundant and the frozen-in condition, 
which inhibits the cross-field plasma flow can perhaps be 
violated. Possible across-the-field flows are the Evershed 
flows in sunspots. Also, various granulation motions can 
perhaps be transverse to the field and hence produce the 
current. Although, in the chromospheric footpoints there 
are only very small random motions detected in horizon¬ 


tal direction (and thus perpendicular to the vertical mag¬ 
netic field), of the order of 10 km s“^, which yields cur¬ 
rents four orders of magnitude below the coronal heating 
requirement. But, on positive side density in the chromo¬ 
sphere is much higher, hence the value of the generated 
current could be enhanced by orders of magnitude. Also, 
further out in the corotating interaction region structures 
llGoslinglll99fi^ . in the solar wind, the magnetic field has a 
transverse to the plasma flow component and hence cur¬ 
rents can be generated. 

We also would like to estimate how ’’efficient” is the 
Sun as a MHD generator: The dimensionless ratio which 
quantifies this, is N = aB'^d/{pUc^) (in CGS units). For 
the most ’’optimistic” estimate (Eq.(9)) along with the 
classical conductivity of cr = 6 x 10^® s“^ for a T = 2 MK 
corona, N = 4.6 x 10® ^ 1. This means that, on one hand, 
the Sun appears to be quite efficient MHD ge nerator, as 

a pical ionised gas MHD generators « 10 I Hevwoo'^ 

. On the other hand, however, despite of this, coronal 
heating budget is still not met by our mechanism. 

Interestingly, our conjecture with MHD generator can 
perhaps also lead to yet another interesting possibility: 
if instead of generating current across the loop (not along 
the loop!) via the flow of plasma across the magnetic field, 
an electric field (created e.g. by some reconnection type 
event) can be applied across ^-coordinate into the plane of 
Fig. 1 i.e. from one edge of the footpoint to the other edge 
of the same footpoint, involving entire bulk of the slab 
(that mimics the coronal loop), t hen the Lorentz force 
reverses its direction (see Fig. 3 in iHevwoo and 

the system can act as a MHD pump, i.e. accelerate the 
plasma flow. Whether this yet another conjecture is of 
any relevance to the solar wind acceleration remains yet 
to be seen. 
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